function dr = ra_to_dr(R, a, B)
dr = a*(2*R+B)/2;
end